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Abstract 

• We prove that the field equations of supergravity for purely time-dependent 

backgrounds, which reduce to those of a one-dimensional sigma model, admit a Lax 
pair representation and are fully integrable. In the case where the effective sigma 
model is on a maximally split non-compact coset U/H (maximal supergravity or 
subsectors of lower supersymmetry supergravities) we are also able to construct a 
completely explicit analytic integration algorithm, adapting a method introduced 
by Kodama et al in a recent paper. The properties of the general integral are 
particularly suggestive. Initial data are represented by a pair Cqi^o where Cq is 
in the CSA of the Lie algebra of U and /iq £ H/VV, is in the compact subgroup 
H modded by the Weyl group of U. At asymptotically early and asymptotically 
late times the Lax operator is always in the Cartan subalgebra and due to the iso- 
spectral property the two limits differ only by the action of some element of the 
Weyl group. Hence the entire cosmic evolution can be seen as a billiard scattering 
with quantized angles defined by the Weyl group. The solution algorithm realizes 
a map from H/W into W. 
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by the Itahan Ministry of University (MIUR) under contract PRIN 2003023852 



1 Introduction 



String Cosmology The cosmological implications of superstring theory have been un- 
der attentive consideration in the last few years from various viewpoints |T]. This in- 
volves the classification and the study of possible time-evolving string backgrounds which 
amounts to the construction, classification and analysis of supergravity solutions depend- 
ing only on time or, more generally, on a low number of coordinates including time. 

The idea of cosmic billiards In this context a quite challenging and potentially highly 
relevant phenomenon for the overall interpretation of extra-dimensions and string dynam- 
ics is provided by the so named cosmic billiard phenomenon [2j,|ni,|l], [Sj- This is based 
on the relation between the cosmological scale factors and the duality groups U of string 
theory. The group U appears as isometry group of the scalar manifold Ai scalar emerging 
in compactifications of 10-dimensional supergravity to lower dimensions D < 10 and de- 
pends both on the geometry of the compact dimensions and on the number of preserved 
supersymmetries Nq < 32. For Nq > 8 the scalar manifold is always a homogeneous 
space U/H. The cosmological scale factors ai{t) associated with the various dimensions 
of supergravity are interpreted as exponentials of those scalar fields hi{t) which lie in the 
Cartan subalgebra of U, while the other scalar fields in U/H correspond to positive roots 
a > of the Lie algebra U. The cosmological evolution is described by a fictitious ball 
that moves in the CSA of U and occasionally bounces on the hyperplanes orthogonal to 
the various roots: the billiard walls. Such bounces represent inversions in the time evo- 
lution of scale factors. Such a scenario was introduced by Damour, Henneaux, Julia and 
Nicolai in a series of papers that we already quoted 05 0, which generalize clas- 

sical results obtained in the context of pure General Relativity 0. In this approach the 
cosmic billiard phenomenon is analyzed as an asymptotic regime in the neighborhood of 
space-like singularities and the billiard walls are seen as delta function potentials provided 
by the various p-forms of supergravity localized at sharp instants of time. 

The smooth cosmic billiard programme In a series of papers [ZI , [HI El lEl involving 
both the present authors and other collaborators it was started and developed what can 
be described as the smooth cosmic billiard programme. This amounts to the study of the 
billiard features, namely inversions in the evolution of cosmological scale factors within the 
framework of exact analytic solutions of supergravity rather than in asymptotic regimes. 
Crucial starting point in this programme was the observation [T that the fundamental 
mathematical setup underlying the appearance of the billiard phenomenon is the so named 
Solvable Lie algebra parametrization of supergravity scalar manifolds, pioneered in JT] 
and later applied to the solution of a large variety of superstring/supergravity problems 
J2] 5 CI] (for a comprehensive review see fSj). Thanks to the solvable parametriza- 
tion, one can establish a precise algorithm to implement the following programme: 

1. Reduce the original supergravity in higher dimensions D > 4 (for instance D = 
10, 11) to a gravity-coupled a-model in D < 3 where gravity is non-dynamical 
and can be eliminated. The target manifold is the non compact coset U/H = 
exp [Solv (U/H)] metrically equivalent to a solvable group manifold. 



1 



2. Utilize the structure of the algebra Solv (U/H) in order to integrate analytically the 
cr~model equations. 



3. Dimensionally oxide the solutions obtained in this way to exact time dependent 
solutions of D > 4 supergravity. 

Within this approach it was proved in [7] that the cosmic billiard phenomenon is indeed 
a general feature of exact time dependent solutions of supergravity and has smooth real- 
izations. Calling h(t) the r-component vector of Cartan fields (where r is the rank of U) 
and ha(t) = a ■ h(t) its projection along any positive root a, a bounce occurs at those 
instant of times tj such that: 



namely when the Cartan field in the direction of some root a inverts its behaviour. All 
higher dimensional bosonic fields (off-diagonal components of the metric g^j_i, or p-forms 
A^'^) are, via the solvable parametrization of U/H, in one-to-one correspondence with 
roots 0Q, <(=^ a, and each bounce on a hyperplane orthogonal to a root a is caused by the 
sudden growing of the field (pa which, in exact smooth solutions, has a typical bell-shaped 
behaviour around the bouncing time t = tt,. 

1.1 Algebraic structure of U/H and integration with the com- 
pensator method 

As alluded to in the above discussion, the actual construction of exact analytic solutions 
performed in j?^ and was based on two ingredients, namely the systematics of dimen- 
sional reduction/oxidation which ultimately allows to reduce supergravity field equations 
to those of a one-dimensional sigma model on a coset U/H and secondly the algebraic 
structure of U/H which allows to device effective integration methods. The essential point 
relative to the second ingredient is that in all supergravity models U/H is non-compact 
and for maximal supersymmetry it is also maximally split, namely the Lie algebra U of U 
is the maximally non compact real section of its own complexification and H is generated 
by the maximal compact subalgebra H C U. In this case the generators ta of EI are 
in one-to-one correspondence with the positive roots a of U and the detailed algebraic 
structure of this geometry was utilized in [7] to device an integration algorithm named the 
compensator method. This latter essentially consists of substituting the original cx-model 
field equations with an equivalent set of differential equations which are more manageable 
since they can be put into triangular form and integrated one-by-one. The geometric 
meaning of these new equations is that they correspond to the conditions to be satis- 
fied by the 6a{t) parameters of an H-gauge transformation L i-^ hh{6) in order for the 
so called solvable gauge of the U/H coset representative L(t) to be preserved, solvable 
gauge being the possibility of writing L(t) as an exponential of the solvable Lie algebra 
S'o/f (U/H). The limitations of the approach developed in [7| |H] and of the compensator 
method which nonetheless was very valuable in producing a rich set of explicit analytic 
solutions are threefold: 



3a G A+ \ 



h„(t) 
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a In its original form, the scheme of [3 |H] apphes only to finite algebras, namely to the 
dimensional reduction of supergravity to = 3, which is the first where all (Bose) 
degrees of freedom can be represented by scalar fields. This is not yet sufficient to 
reveal the full algebraic structure underlying time evolution and billiard dynamics. 
Indeed to discuss this latter one can still dimensionally reduce to D = 2 and D = 1 
where the algebra U is extended to an affine or hyperbolic Kac-Moody algebra. It 
is therefore mandatory to consider the field theoretical mechanisms of such Kac- 
Moody extensions and subsequently extend the so far considered algebraically based 
integration methods to Kac-Moody algebras. 

b Heavily relying on the maximally split character of the pair {U , H} the original form of 
the scheme presented in j3IH] is not directly applicable to the case of lower supersym- 
metry Nq < 32 and hence to more realistic and more interesting compactifications 
than pure toroidal ones. 

c The compensator method, although very valuable in obtaining explicit solutions, re- 
places one differential system with another one and therefore does not answer the 
question whether the system is a fully integrable system nor it is able to provide 
a uniformly given, general integral, depending on as many integration constants 
as necessary. Indeed the solution strategy is based on triangulization of the set of 
equations and the choice of parametrization of the H-subgroup element h (order of 
the rotations exp [9ait)ta] ) in such a way that the resulting differential system for 
the rotation angles 9ait) be triangular is something that has to be decided case by 
case and both by intuition and by trial and error. 

The two recent papers [10^ and [9J have presented some substantial advances with respect 
to points a] and b] of the above list. 

a] : In [lOj it was clarified the general field theoretical mechanism underlying Kac- 
Moody extensions in supergravity and the general pattern of such extensions in the re- 
duction to D = 2 of all Z) = 4 supergravity theories with all number of supercharges 
Nq was also presented. This means that such extensions are now under control for all 
compactifications of both M-theory and string theory on all type of compact geometries. 
Indeed these latter lead to the various different types of -D = 4 supergravities. Although 
integration methods for the Kac-Moody analogue of the cr-model equations were not pre- 
sented in [TU], yet the results there obtained provide the necessary prerequisite in the 
development of such a programme. Indeed the field theoretical mechanism of the exten- 
sion is what one needs to consider in order to rewrite the actual D = 2 or D = 1 field 
equations in a manifestly Kac-Moody covariant way. 

b] : In [S], instead, exploiting in a systematic way the Tits Satake theory of real forms 
[TB] for simple algebras and the Tits-Satake projection of root-spaces, a decisive advance 
was made with respect to point b] of the above list. Indeed it was shown how the 
differential equations of the a-model on a non maximally split non compact coset U/H, 
typically emerging from non maximal supergravities, can be consistently reduced to those 
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of a (T-model on a maximally split one: Uts/Hts- By Vts C U we denote the Tits- 
Satake subalgebra of the original algebra U which is uniquely defined by the involutive 
automorphism a associated with the real form of U. By Mts C Uts we denote the 
maximal compact subalgebra of the Tits Satake subalgebra which, as in all maximally 
split cases, has a number of generators equal to the number of positive roots in the Au 
root system of Uts- This latter is the Tits Satake projection of the U-root system: 

n : Au Au (1.2) 

obtained by setting to zero all transverse components a± of each root a. Indeed the 
automorphism o"^ = 1 defining the real form of U splits the Cartan subalgebra CSA C U 
in two eigenspaces Ti.""'^' and 7Y'^°"*p, respectively corresponding to the eigenvalues ±1. 
Each root a, which is a linear functional on the CSA, is accordingly decomposed into a 
parallel and into a transverse part to the non compact Cartan subalgebra Ti'^ '^ -. 

a = aii®a± ; (^"-'O = ; a,, ^ 0. (1.3) 

The Tits Satake projection projects each root a onto its parallel part a\\. 

As already mentioned it was shown in |U] that in the non maximally split cases one can 
consistently reduce the cr-model equations of U/H to those of Uts/Hts and furthermore 
it was proved that one can lift back any solution of these latter to a full fledged solution of 
the original system by means of the action of an automorphism compact group named the 
paint group Gpamt C U with respect to which the entire solvable algebra Solvu transforms 
as a representation: 

[Gpamt, SolVu] = SolVu. (1.4) 

Hence a large class of solutions of the cr-model on non-maximally split coset manifolds, 
and presumably the essentially relevant class for billiard dynamics, since it contains all 
the Cartan fields and all the walls on which the cosmic ball can bounce, can be obtained 
from the general solution of a-models on maximally split cosets. 

1.2 The present paper solves point c]: full integrability via LAX 
representation 

In the present paper we shall present a complete solution for the third point in the 
above list of problems and required generalizations. Indeed we shall prove that the one- 
dimensional cr-model on maximally split cosets is completely integrable in the sense of 
integrable systems, since the first order equations for the tangent vector to a geodesic can 
be recast in the classical form of a LAX system. Furthermore, adapting to our case an 
algorithm developed in il7|, we can write a closed analytic form for the general integral, 
depending on a complete set of integration constants. 

This new algorithm replaces the compensator method and does not require the solution 
of neither differential nor algebraic equations. The solution of the equations for the time 
parameter ranging from — oo to +00 is directly constructed in terms of the initial data 
fixed say at t = 0. It is particularly interesting to study the general properties of the 
solution algorithm. As we shall explain, the initial data parametrizing the value Lq of the 
Lax operator L{t) at t = can be represented as a pair: 
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1. An element of the Cartan subalgebra Co G CSA C U. 

2. An element ho eB. = exp [H] of the maximally compact subgroup H C U. 
The solution algorithm Sk provides a general integral 

Sk ■■ {Co, ho} =^ L(t\Co,ho) (1.5) 
with the following properties: 

1. V w G W where W C H C U is the Weyl group of the Lie algebra U 

L (t\Co,w ■ ho) = L (t\ ■ Co ■ w, ho) . (1.6) 

2. The limits of the Lax operator L {t \ Co, ho) aX t = ±oo lie in the Cartan subalgebra: 

lim L{t\Co,ho) = L±^{Co,ho) G CSA. (1.7) 

3. At any instant of time the eigenvalues {Ai A2 . . . A„} of the Lax operator are the 
same and are given by 

A. = w, (Co) (1.8) 

where Wj are the weights of the representation of U to which the Lax operator is 
assigned. 

Property (2) and property (3) combined together imply that the two asymptotic values 
L-too {Co, ho) of the Lax operator are necessarily related to each other by some element 
i7 G W of the Weyl group which represents a sort of topological charge of the solution: 

3 G W : {Co, ho) = ■ (Co, ho) -n. (1.9) 

Indeed since the Weyl group is discrete, by varying the initial datum ho a solution which is 
characterized by a charge Qi, cannot continuously be deformed into another characterized 
by a different charge ^2- 

It happens furthermore that the topological charge of a solution is independent of 
the choice of Cq G CSA and just depends on ho- We conclude that the solution algorithm 
Sk developed in [T7j induces a map: 

^ ^ ^ >V = ^1 (1-10) 



from the non-simply connected coset manifold ^ to its own homotopy group. 

Let us now proceed to the derivation of these results and to the illustration of their 
application to cosmic billiards. 
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2 One-dimensional sigma models on maximally split 
coset manifolds and Lax pairs 

In view of the discussion presented in the introduction our mathematical object of study 
is reduced to be the one- dimensional a-model on a maximally split hermitian symmetric 
coset manifold U/H. The corresponding action principle is given by: 



where hjj is the canonical metric on the coset, which we parameterize by scalar fields 
named 0^. By definition, the coset U/H is hermitian symmetric if the Lie algebra U of 
the group U and the subalgebra HI of the subgroup H fulfill the relations: 

[M, H] C H, [H, K] C K, [K, K]cM (2.2) 

having called IK the orthogonal complement of HI in U, namely: 

HcU, U = H©K. (2.3) 

In order to asset the further property of being maximally split we assume that U is a simple 
Lie algebra in its maximally non-compact real section and we normalize the Cartan-Weyl 
commutation relation in the standard form: 

[Hi , Hj] = 0, 

[E" , E-"] = a ■ n, 

[E"^,E^^] = Ar„,,,^, (2.4) 

where Hi denotes a basis for the Cartan subalgebra (CSA), the root vectors of the 
root system A and i?"' the corresponding step operators. Then U/H is maximally split if 
the subalgebra HI is chosen to be the maximally compact subalgebra which has dimension 
equal to the number of positive roots and which can be defined as: 

H = Span{t^} = Span{(E" - E''')} . (2.5) 

Eq. (I2.5|l defines the subalgebra and also a basis for it given by the generators = 
{E" — E~"). In this case the complementary space IK is given by 

K = Span{K^} = Span{Hi, ^(S" + E'")} (2.6) 

v2 

which, similarly to eq. (|2.5|) introduces a canonical basis of generators Ky^. 

For any parametrization of the coset provided by a coset representative L(0) one can 
write down the U-algebra valued left invariant one-form 

= L"^ = + = W + V (2.7) 
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where V is the coset manifold vielbein using which one can rewrite action ()2.1|) as follows 

S= [ dtTi (VV) . (2.8) 



In eq. ()2.8j) it is understood that we have chosen a linear representation of the Lie 
algebra U (which one is irrelevant, typically the lowest dimensional one) and therefore 
the vielbein has become a n x n matrix valued one-form, n being the dimension of the 
chosen representation: 

V{t) = Ka V^j (0(t)) ^ = . (2.9) 

Our goal is to solve the classical field equation of the model ()2.ip . namely calculate the 
geodesies of the manifold U/H. Let us calculate the variation of the action 

5S = 2 j dt TiiVSV). (2.10) 

In order to derive a convenient expression for the variation of the vielbein SV which enters 
eq. ()2.10|) . we represent the variation of the left-invariant one-form Q defined in eq. ()2.7p 
in the following equivalent way: 

6n = [n, h-^h] + ^(L-i(5L) (2.11) 

and take into account that 

L-^(5L = 5cu"t^ + 6v^Ka = 6u + 6v, (2. 12) 

where 6v^ and 6u°' are functionals of the variation of the fields 6(j)^ , and there is the 
one-to-one correspondence between 6v^ and 6(f)^ (and between and (j)^ modulo trivial 
constants), so that the transition from 6(j)^ to 6v^ can be treated as a change of the basis. 
Substituting the relations (|2.7|) and (|2.12|) into eq. (|2.11|) we obtain the following new 
relation: 

SW + SV=[W + V, Su + Sv] + j^{5uj + 5v). (2.13) 

Now, using eq. ()2.2|l it is easy to extract from eq. ()2.13j) the projection along the coset 
generators and derive 5V 

dV = [W, 6v] + [V, 6uj] + ^6v. (2.14) 

Substituting this expression into ()2.1()|1 one can easily obtain the variation of the action 
5S 



5S = 2 j dt Ti (v[W, 5v] + V[V, 5uj] + V^Sv 

= 2 J dtTr (^(^[V, W] - j^V^ 6v^ . (2.15) 



The equations of motion follow immediately from the vanishing of the action variation 
5S = 



-V=[V,W] (2.16) 

in the form of a Lax pair representation. Hence the resulting system is indeed integrable. 
The direct consequence of the Lax pair representation ()2.16|) is the existence of the fol- 
lowing integrals of motion: 

Ik = Tr {V'') . (2.17) 

for all integers k < n, having called n the dimension of the matrix V namely the dimension 
of the chosen linear representation of the Lie algebra U, as we have already specified. 



2.1 Nomizu connection for maximally split cosets and the Lax 
pair representation 

Now in order to make contact with our previous work and with the discussion of cosmic 
billiards, we want to compare the Lax pair representation of the field equations (j2.16j) 
with the form of the same equations derived in |7j by means of the Nomizu connection on 
solvable Lie algebras. There showed that the field equations of the purely time dependent 
(j-model ()2.1|) . which is what we are supposed to solve in our quest for time dependent 
solutions of supergravity, can be written as follows: 

+ ricY^Y^' = (2.18) 

where Y^ denotes the purely time dependent tangent vectors to the geodesic in an an- 
holonomic basis: 

A _ \Y' = VI {(f)) i e CSA 



\ Y" = \/2 Vf {(j)) (j)^ a G positive root system A+ 

Vf" {(f>) d(j)^ being the vielbein of the target manifold we are considering, namely the same 
object already introduced in the previous section. In eq. ()2.18|1 the symbol T^(^ denotes 
the components of the Levi-Civita connection in the chosen anholonomic basis. Explic- 
itly they are related to the components of the Levi-Civita connection in an arbitrary 
holonomic basis by: 

r^c = ^mV^ViV^ - dK{Vf)ViV^ (2.20) 
where the inverse vielbein is defined in the usual way: 

y/yi = 5^ (2.21) 

The basic idea of which was exploited together with the compensator method in order 
to construct explicit solutions, is the following. The connection Vj^fj can be identified 
with the Nomizu connection defined on a solvable Lie algebra, if the coset representative 
L from which we construct the vielbein is solvable, namely if it is represented as the 
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exponential of the associated solvable Lie algebra So/f (U/H). In fact, as we can read in 
[TH] once we have defined over Solv a non degenerate, positive definite symmetric form: 

( , ) : Solv (g) Solv — > M, 
{X,Y) = {Y,X) (2.22) 

whose lifting to the manifold produces the metric, the covariant derivative is defined 
through the Nomizu operator: 

VX e Solv : ILx : Solv — > Solv (2.23) 

so that 

VX, Y,Z e Solv : 2{Z , hxY) = {Z, [X, Y]) - {X, [Y, Z\) - (F, [X, Z\) (2.24) 
while the Riemann curvature 2-form is given by the commutator of two Nomizu operators: 

i?^^(X,F) = (ly, {[Lx,Ly] - L[x,y]}^) • (2.25) 
This implies that the covariant derivative explicitly reads: 

LxF = rfyZ (2.26) 

where 

Yly = i ((z, [X, Y\) - (X, [y, z\) - (y, [x, z\)) vx, y, z g Soiv . 

(2.27) 

Eq. ()2.27|) is true for any solvable Lie algebra, but in the case of maximally non- 
compact, split algebras we can write a general form for T^yy namely: 



= 0, 



a/3) 



J- a/3 - Pa - 2 

^a+pp = ^pa+P = ^^ap (2.28) 

where X"^ is defined by the commutator [E^ , Ep] = Nap E^+p, as usual. The explicit 
form ()2.28|) follows from the choice of the non degenerate metric: 

{Hi , Hj) = 2 Sij, 
{Hi , Ea) = 0, 

{Ea^Efs) = 6a,p (2.29) 
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V7-^j, Hj G CSA and V-Eq,, step operator associated with a positive root a G A+. 

Hence in the case of maximally split algebras the first order equations, take the general 
form: 



+ N^pY'^y'"^'^ - c^i Y'Y'" = (2.30) 

/3eA+ 

which follows from eq. ()2.28|1 . 

Our next point is that equations ()2.3()j) are identical to the Lax pair equations ()2.16j) 
if we identify: 

V = v^Ka = vni + y^liE'^ + E-"), 

On one hand insertion of (|2.31|) into (j2.16|) . upon use of the standard Cartan Weyl commu- 
tation relations ()2.4|) . yields ()2.30p . on the other hand eqs. ()2.31|) are justified as follows. 
From eq. ()2.19|) we have 

l^" = ^y" (2.32) 
V2 

while, on the other hand we have that the condition for the coset representative h{(f)) to 
be solvable is given by: 

V'' = V2W''. (2.33) 

Together eqs. ()2.32|) and ()2.33|) imply ()2.3H) . In this way we have proved that the first 
order equations ()2.30|) studied in jTHHUH] and for which we have also recently found bounc- 
ing solutions also in the interesting case of F4(4)/Usp(6) x SU(2), which is relevant for 
= 6 supergravity, are actually fully integrable since they admit the Lax pair rewriting 
()2.1(ij) . In [3ISlin] eqs. ()2.30j) were addressed by means of the compensator method. Let 
us recall its form in order to compare it with the general integral formula provided by the 
algorithm of p!7j . 



2.2 Minireview of the compensator method 

As we have seen above, the condition which ensure that a coset representative L(0) is 
solvable, namely that can be represented as the exponential of the solvable Lie algebra 
Solv{\J /H) is given by the linear relation ()2.33p on the components of the left-invariant 
one-form ()2.7|1 . The compensator method streams from the following argument. Given 
a solvable coset representative L(0), one asks the question how many other solvable 
representatives are there in the same equivalence class. This amounts to investigating 
the condition to be satisfied by the H-gauge transformation 

L ^ L/i = L, 

h = exp (2.34) 
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in order for the solvable gauge ()2.33|) to be preserved. This condition is a system of 
differential equations, namely: 



x/2 



tr {h-\e) dh{e) Q = V" (-A(^)/ + /^(^)/) + V (2.35) 



In eq. (j2.35p the matrix A{9) is the adjoint representation of /i G H and D{9) is the 
D-representation of the same group element which acts on the complementary space K 
and which depends case to case 

h-^KAh = D{e)/^B ■ (2.36) 

The differential system (j2.35|) is actually equivalent to the original system (j2.3(Jj) . To see 
this it suffices to argue as follows. A simple solution of the first order equations (j2.3U|) 
is easily obtained by setting Y"" = and = d = const, namely we can begin with a 
constant vector in the direction of the CSA. Such a solution is named the the normal form 
of the tangent vector. In the language of billiard dynamics it corresponds to a fictitious 
ball that moves on a straight line with a constant velocity. All other solutions of eqs. 
(|2.3U|) can be obtained from the normal form solution by means of successive rotations of 
the compact group, with parameters 9[t] satisfying the differential equation ()2.35|) . The 
advantage of this method, emphasized in [7] where we introduced it is that we can choose 
rotations in such a way that at each successive rotation we obtain an equation which is 
fully integrable in terms of the integral of the previous ones. How to do this, however, can 
not be prescribed in general and furthermore it is not obvious how many steps one can do 
of this type. In other words, although the compensator method provides a valuable tool 
to obtain several explicit solutions, it does not suffice to reveal the complete integrability 
which is instead revealed by the Lax pair representation ()2.16p . 



2.3 The A2 case with the compensator method 

To appreciate the bearing of the compensator method relative to the general integration 
algorithm that we shall present in next section let us consider the simplest example of 
non abelian maximally split coset, already studied and solved in [7j. This is the manifold 

U/H = SL(3,M)/SO(3) (2.37) 

and the Lie algebra U = SL(3,M) is the maximally non-compact real form of the simple 
complex Lie algebra A2. The coset decomposition (j2.3p takes the form: 

SL(3,R) = S0(3) ©K (2.38) 

where 

SL(3, M) = SpaniHi, E", E""}, i = 1, 2, a = 1, 2, 3, 

K = Span {Ka} = Span S^Hu {E" + E^") | , 

S0(3) = Span {t^} = Span { {E" - E'") ] . (2.39) 
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In the fundamental 3x3 representation the generators of SL(3,]R) are given by: 



1 
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-2/ 




/o 





l\ 




















0/ 




/o 





0^ 













. (2.40) 












In [7j we considered the differential system ()2.35p generated by taking as generating 
solution in normal form the following: 



1^1 



(z/i,z/2, 0,0,0), 

UJ — K 3u! + K 

^2 



4^/2 



4^6 



(2.41) 



where u and k, are just two alternative constant parameters replacing the components 1^12 
of the normal solution along the two Cartan generators. Parametrizing the h £ SO (3) 
gauge group element as follows: 

h = exp [03 (i) ts] exp [O^it) t^] exp [9,{t) t^] (2.42) 

we found in |7| that the differential system ()2.35|) takes the form 

= Ut) -\ujsm{2 9s{t)), 

= se2{t)-[K + uj cos(2 e^it))] sin(2 62(1)), 

= 169i{t) + [k + K cos{2 92{t)) + uj[cos{202it)) - 3]cos{29^{t))] sin(2 0i(t))- 

-8ujsm^{9i{t)) sin(02(t)) sin(2 03(t)). (2.43) 

This system has the requested triangular property and the first two equations are imme- 
diately solved by 



arccos 



arctan 




(2.44) 



where A2,3 are two integration constants. The third equation of ()2.43|) was also solved 
in [7] with some ad hoc elaborations and in this way we obtained a general integral 
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depending on the five integration constants {uj, Ai, A2, A3}. This is obviously consistent 
with the full integrability of the system demonstrated by the Lax pair representation, 
but still required some ad hoc elaborations. Furthermore for higher rank algebras like for 
instance ^3, various experiments showed that although the system ()2.H5j) can be arranged 
to triangular for a few rotations, it is not possible to make it completely triangular if we 
include rotations along all the generators of H. These limitations are now completely 
overcome by the general integration formulae which we present in the next section. 



3 The general integration algorithm. 

For maximally split cosets U/H, choosing a solvable coset representative L(0), there is 
an established algorithm jT^j to construct the general solution of eq. ()2.1(i|l . For the 
maximally non compact real section of any simple Lie algebra U there always exists a 
permutation matrix O such that in the Weyl basis the Lax pair representation ()2.16|) can 
be brought to the following general form 

ft^ = [L,P] (3.1) 

where 

L = OVO^, P = OWO^ (3.2) 

and the matrix P is given by the following projection of L 

P = n(L) := L>o - L^o, (3.3) 

L>o (<o) denoting the strictly upper (lower) triangular part of the n x n matrix L. The 
existence of the matrix O, such that eq. ()3.3|) holds true is guaranteed by a theorem 
proved in Helgason's book jTHI, by the maximally split nature of the coset U/H and 
by the condition of solvability ()2.33p imposed on the coset representative. The quoted 
theorem states that for any linear representation of a solvable algebra there always exist 
a basis where all generators are upper triangular matrices. Invoking this theorem for 
the solvable algebra SolviXJ/B.) generated by {Hi , -E"} we conclude that we can always 
bring the semisimple operators Hi to diagonal form and the nilpotent operators to 
strictly upper triangular form. In the same basis the step-down operators which are 
the adjoints of the corresponding E'^ are necessarily strictly lower triangular. Since the 
connection W is by definition given hj W = ta W'^ = W" (i?" — E~"), while V is given 
hj V = Hi + V" {E"' — E~"), we see that in the basis where E°' is upper triangular 
and Hi diagonal, eq. (|3.3|) is nothing else but the solvability condition (j2.33p . 

This established, we can proceed to apply the integration algorithm jT7j in order 
to write the general integral of the Lax equation ()3.1|) . Actually this is nothing else 
than an instance of the inverse scattering method. Indeed equation p.ip represents the 
compatibility condition for the following linear system exhibiting the iso-spectral property 
of L 

^ = (3.4) 

dt ^ ' 
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where is the eigenmatrix, namely the matrix whose i-th row is the eigenvector (f{t, Aj) 
corresponding to the eigenvalue Aj of the Lax operator L{t) at time t and A is the diagonal 
matrix of eigenvalues, which are constant throughout the whole time flow: 

= [V^(Ai),...,^(A„)r=[^,(A.)]^<^_^.<^, 
A = diag(Ai,...,A„). (3.5) 

The solution of ()3.4|1 for the Lax operator is given by the following explicit form of the 
matrix elements: 



(3.6) 



k=l 



The eigenvectors of the Lax operator at each instant of time, which define the eigenmatrix 
^(i), and the columns of the its inverse \l/~^(t), are expressed in closed form in terms of 
the initial data at some conventional instant of time, say at t = 0. 
Explicitly we have: 



ipi{Xj,t) 



1pj{\i,t) 



-Ait 



y/Dj{t)Dj_,it) 



Det 



Cll 



-1 ^?(A,)A 



Det 



Cll 



\ 



^°(A.) ) 



(3.7) 



where the time dependent matrix Cij{t) is defined below: 



N 



,.(t) = J]e-2^'=V°(Afc)^J(A.) 



(3.^ 



and 



k=l 

^°(Afc) 



V5j(Afc,0) 
^i(Afc,0) 



(3.9) 



are the eigenvectors and their adjoints calculated at t = 0. These constant vectors consti- 
tute the initial data of the problem and provide the integration constants. Finally Dkit) 
denotes the determinant of the k x k matrix with entries Cij{t) 



Dk{t) = Det 



l<i,j<k 



(3.10) 



Note that Cjj(O) = 6ij and Dk{0) = 1 as well as -Do(^) := 1- 
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3.1 Parametrization of the initial data 



Let us now reconsider the integration formulae presented in the previous subsection keep- 
ing in mind that, at any instant of time t, the Lax operator L{t) is an element of the Lie 
algebra U lying in the orthogonal complement K to the subalgebra H. This Lie algebra el- 
ement is expressed in the chosen n-dimensional representation. Diagonalizing L{t) means 
nothing else than bringing it inside the Cartan subalgebra CSA C K C U, which can 
always be performed by conjugation with an element of the maximal compact subgroup. 
Hence at every instant of time there exists an element h{t) G exp [H] such that: 

3 h{t) G exp [M] \ h-\t) L(t) hit) = Co G CSA. (3.11) 

The n X n matrix Co is by construction diagonal and has therefore on its diagonal the 
above mentioned constant eigenvalues A^. What are these? They are just the n-weights 
of the chosen n-dimensional representation evaluated on the Cartan element Cq. Calling 
Wj the weights, we just have: 

\i = Wi (Co) =Wi-l/ {i = l,...,n) (3.12) 

having denoted by the r-vector 'P' the r components of the Cartan element Co in a basis 
of the Cartan subalgebra: 

Co = UiW 2 = 1, ...,r = rank(U). (3.13) 

Equation p.llj) is in particular true at the initial time t = 0, where there exists some 
constant ho G exp [H] such that 

/lo ' ^(0) ho = Co = diag (Ai . . . , A„) . (3. 14) 

At this point all items entering the solution algorithm described in the previous section 
obtain their proper group-theoretical interpretation. The diagonal eigenvalue matrix A is 
nothing else but the CSA element Co evaluated in the chosen n-dimensional representa- 
tion. The eigenmatrix ^1^(0) which provides the initial data of the problem is nothing else 
but the if-group element ho. This explains what we anticipated in section IT!^ when we 
stated that the initial data are a pair of CSA algebra element and a H-subgroup element. 
We still have to verify the other properties relative to the Weyl group. These are also 
easily understood. By definition the Weyl group W is that discrete subgroup of H which 
maps the Cartan subalgebra into itself: 

Vw G W, VC G CSA w-^-C-w G CSA. (3.15) 

Hence if we choose ho inside the Weyl group the initial value of the Lax operator Lo = L{0) 
remains a diagonal matrix in the Cartan subalgebra. If Lo is diagonal the connection 
P(0) = vanishes at the origin and from Lax equation we obtain 

^L\t=o = (3.16) 
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which imphes that the solution is constant throughout the whole flow: L{t) = L{0). 
This suffices to prove eq. ()1.6|) . On the other the so called sorting property, namely the 
statement ()1.7j) was proved by Kodama et al in As already anticipated in section lL2l 
eq. p.fij) . together with the iso-spectral property, namely with the form of the eigenvalues 
()H.12j) implies that eq. has also to be true. Finally the independence of the charge 

Q from the eigenvalues, namely from the choice of the CSA element Cq, together with 
eq. p.6|) implies that the map induced by the solution algorithm is indeed of the type 
described in eq. p.lO|) . 



4 The As example with numerical plots 

In order to exemplify the value and the use of the general integration algorithm, we will 
apply it to the case of the manifold 



M:- 



SL(4,B 



(4.1) 



which is maximally split of rank r = 3 and corresponds to the Lie algebra ^3. As in the 
previous A2 case we can use the fundamental representation as working representation for 
the Lax operator, which means that we deal with 4x4 matrices. The dimension of the 
coset is 9 and there are 6 positive roots: 



"5 



{72,0,0} 





1 

V2' 



2 _2_ 

3' 

1 2 



^2 



L _ /3 

V2' A/'J'^ 



1 

1 

V2' 



^ 

2' " 



1 2 



(4.2) 



The matrices representing the Cartan generators well adapted to the representation (j4.2j) 
are the following ones: 



Hi 











o\ 





1 

V2 




















V 








0/ 



Ho 



I 2V3 




\ 





1 

Ve 








o\ 







1 
















[2 

V 3 





V 











0/ 





1 

2^3 










1 

2 





\ 





-V3 



while the 6 step-up operators associated with the positive roots are given by: 



E 



/O 1 0\ 





\0 0/ 



E 



a2 



/O 0\ 

10 



\o oy 



(4.3) 



(4.4) 



16 



/o 








0\ 




/o 





1 


0\ 





































1 






















V 












0/ 


fo 








0\ 























1 










































^0 








V 












0/ 



(4.5) 



(4.6) 



The step-down operators associated with the negative roots are just the transposed of the 
corresponding step-up operators. 

The H subgroup is 0(4), namely the set of all orthogonal 4 x 4-matrices. The Weyl 
group W C 0(4) C SL(4, M) is S4 and its elements are represented by the orthogonal 
matrices associated with the 4! permutations P e 5*4 in the following way: 



(Op),P(,) 

So for instance the cyclic permutation 



= 1, 

= otherwise. 



(4.7) 



P = 



12 3 4 

J, J, J, J, 

4 12 3 



(4.8) 



is represented by the matrix 



/ 1 \ 

10 

10 

\ 1 / 



(4.9) 



Having defined the model we can now study the solutions produced by the general inte- 
gration algorithm implementing it on a computer and choosing numerical values for the 
initial data. 

Let us begin with the choice of the Cartan subalgebra element Cq. Setting: 



i/i = - 



Ai + A2 -h 2 A3 
V2 



We obtain: 



Co = ^2 



i=l 





A1 + 3A2 






; = 


\/6 






Ai — A2 


-A3 





\ 





A3 














A2 

















-2Ai 



(4.10) 



(4.11) 
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We are free to choose the eigenvalues Aj. We exhibit the solution for the simple choice: 



A. 



Co 



/-6 0\ 

3 

2 

\ 1/ 



(4.12) 



Next we have to choose the iJ-subgroup element. Our numerical choice is given by the 
following orthogonal matrix: 



8 

-(-1+V3) 

4v/2 
-(3+2v^) 

-8+v^ 
8^2 



-5 (-1+V3j 
16 

V2+3x/6 
16 

7+4 
16^2 

-(lo+Vs) 

16^2 



l+\/3 
4^2 

1+v^ 
4 

-6+V3 
8 



-1-5 
16 

-(-7+V^) 

8v^ 

-4+v^ 
16 

11+2 
16 



(4.13) 



which is just a product of rotations along the six generators of 0(4) with angles that are 
multiples of 7r/3 or 7r/4 in order to obtain simple but non trivial entries all over the place. 
With such initial data the computer programme calculates the Lax operator L{t) at all 
times and finds that the limits at ±oo are as follows: 



lim L{t) 



lim L{t) 



/3 \ 

2 

10 

\0 -6/ 

/-6 0\ 

10 

2 

\ 3/ 



(4.14) 



This means that the image of ho G 0(4) under the map Vk of eq. (jl.lO|) is the permutation 
displayed below: 



K 



0(4) 



3 hn 



TTi {ho 



/O 1\ 

10 

10 

\1 0/ 



(4.15) 



The permutation of the eigenvalues occurs through a smooth process that realizes the 
billiard bounces and that can be appreciated through the inspection of plots of the Cartan 
fields. The most instructive plots are those of the differential and integrated Cartan fields. 
Given the explicit solution for the Lax operator L{t) obtained through implementation of 
the integration algorithm, let us define: 

t.At) = Tr[K.L(t)] (4.16) 
the components of the tangent vector along the geodesic. Let us also define the projection 
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Figure 1: Plots of the differential and integrated Cartan fields Xaiit) and ha^{t), along the 
first simple root of the Lie algebra and for the solution generated by the initial data 
Co and h^, given in eqs. ()4.12|) and ()4.13|) . respectively. The plot on the left is ha^it), 
while the plot on the right is Xaiit). We see the bounce at t = 0.29488. 




Figure 2: Plots of the differential and integrated Cartan fields Xa2(t) and ha2{t), along 
the second simple root of the Lie algebra and for the solution generated by the initial 
data Cq and ho, given in eqs. ()4.12p and ()4.13|) . respectively. The plot on the left is ha2{t), 
while the plot on the right is Xa2(t)- We see the bounce at t = —0.416 
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Figure 3: Plots of the differential and integrated Cartan fields Xasit) and h^^it), along 
the third simple root of the Lie algebra and for the solution generated by the initial 
data Cq and Hq, given in eqs. ()4.12p and ()4.13|) . respectively. The plot on the left is ha^it), 
while the plot on the right is Xasit)- We see the bounce at t = 0.160124. 



of this tangent vector along the simple roots 

XaAt) = l^it)-a, ; 2 = 1,2,3. (4.17) 

Finally the Cartan fields, that in supergravity correspond to the cosmological scale factors 
and to the dilatons ^, are defined by simple integration: 

h^Xt) = [ XaA^) di ; z = 1,2,3 . (4.18) 

J to 

In fig.s lll2l31 we present, for the considered explicit solution the plots of Xoiit) and ha.{t). 
Billiard bounces correspond to maxima or minima of the simple root integrated Cartan 
fields ha.{t) or, equivalently to zeros of the differential ones Xoiit). Inspection of the plots 
reveals that the considered solution, characterized by the topological charge of equation 
()4.15j) admits three bounces, respectively located at: 

t = 0.160124 ; t = -0.416 ; t = 0.29488 . (4.19) 

5 Billiard interpretation 

Having in mind the explicit example we have just illustrated we can now emphasize the 
physical billiard interpretation of the mathematical results we obtained. 

"'^See [3 Ej for detailed explanations about the oxidation mechanism that lifts the tr-model solutions 
to full-fledged solutions of supergravity theory in higher dimensions. See also and [5] for the relation 
between the cr-model and D = A supergravity with all number of supercharges. 
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At those times when the Lax operator L{t) lies in the Cartan subalgebra, the fictitious 
cosmic ball moves on a straight line, the components of L{t) along the CSA generators 
Tix being the components of its velocity: 

If^it) = Tr [n.L{t)] . (5.1) 

This means that the multidimensional Universe has a constant rate (positive or negative) 
exponential expansion (contraction) in the corresponding dimensions. The asymptotic 
theorem states that for finite Lie algebras, namely as long as we confine our attention to 
time-dependent backgrounds that can be retrieved by oxidation from D = 3 supergravity 
(see [ZllHlinilini), the Universe is in such a constant velocity state at asymptotically early 
and asymptotically late times. Furthermore given the initial state: 

1^,-00 > (5.2) 

at t = —00 (Big Bang time) the entire cosmic evolution can be seen as a smooth process 
that brings the Universe to another such state: 

+00 > . (5.3) 

The surprising result is that there is a finite countable numbers of such possible states, 
as many as the elements of the Weyl group W. Indeed the final Universe velocity is 
necessarily the image of its initial one under the action of some element of W. The 
transition from one state to the other occurs through one or several bounces on the 
dynamical walls of the Weyl chamber that as we have already emphasized in previous 
papers raise and decay smoothly. 



6 Conclusions and Perspectives 

We could now show explicit analytic formulae for the general integral of simple models 
like the A2 model. It is very easy to obtain them by means of computer codes written in 
MATHEMATICA but they are very lengthy and clumsy to be displayed so that we prefer 
not to. The relevant point is that the integration algorithm is fully explicit and we hope 
we have been able to illustrate this fact. This enables us to apply our method to any 
case of interest, as long as the target manifold is maximally split. We can also generate, 
by means of the paint group rotations entire classes of solutions in the non maximally 
split case. However they do not yet constitute a general integral in these cases. The open 
problem is therefore that of extending the results presented in this article in the following 
directions: 

1. For the non maximally split non compact cosets appearing in supergravity with non 
maximal SUSY. 

2. For the infinite dimensional Kac-Moody algebras generated by reduction to D = 2 
and D = 1. This is particularly important in view of the results we have achieved 
on asymptotic regimes. It seems quite clear that the chaotic Kasner behaviour 
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predicted around initial or final singularities must be related to the essential change 
in structure of the Weyl group for Kac-Moody algebras. The presence of infinitely 
many roots and moreover the presence of space-like or null-like roots makes the 
Weyl group much more complicated and apt to create smooth solutions with singular 
behaviour if there is an extension of our main results to the Kac-Moody case. 

3. For (T-models deformed by the presence of a potential, like it happens when fluxes are 
turned on and one considers gauged supergravity rather than ungauged supergravity. 

4. Generalization of the results obtained in this paper to the case of one-dimensional 
super-coset sigma models, which comprise both bosonic and fermionic fields. 

On all of these four issues we are presently working and hope to be able to present new 
conclusions soon. 
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